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Abstract 

Effects of the 1-3 leptonic mixing on the solar neutrino observables are studied and the 
signatures of non-zero 613 are identified. For this we have re-derived the formula for 3z^- 
survival probability including all relevant corrections and constructed the iso-contours of 
observables in the sin^ ^12 — sin^ ^13 plane. Analysis of the solar neutrino data gives sin^ ^13 = 
0.007l[J;[!g? (90% C.L.) for Am^ = 8 • 10"^ eV^. The combination of the ratio CC/NC at 
SNO and gallium production rate selects sin^ ^13 = 0.017 it 0.026 (Icr). The global fit of all 
oscillation data leads to zero best value of sin^ ^13. The sensitivity (Icr error) of future solar 
neutrino studies to sin^ ^13 can be improved down to 0.01 - 0.02 by precise measurements of 
the pp-neutrino flux and the CC/NC ratio as well as spectrum distortion at high (E > 4: MeV) 
energies. Combination of experimental results sensitive to the low and high energy parts of 
the solar neutrino spectrum resolves the degeneracy of angles ^13 and 612- Comparison of 
sin^ ^13 as well as sin^ ^12 measured in the solar neutrinos and in the reactor/accelerator 
experiments may reveal new effects which can not be seen otherwise. 

1 Introduction 



Determination of the 1-3 mixing parameterized by the angle 6is is of great importance for phe- 
nomenology, for future experimental programs, and eventually, for understanding the underlying 
physics. 

At present, for sin^ 6'i3 we only have the upper bounds. The CHOOZ [1] reactor experiment 
(see also Palo Verde [2]) gives 

sin^^i3 < 0.036, 90% C.L. (1) 

for Am^3 = 2.5 ■ 10~^ eV^. The recent global fit of the oscillation results which includes also the 
solar salt phase II SNO, KamLAND and atmospheric neutrino data leads to [3] 

sin^ ^13 < 0.026 (0.047), 90% (3a) C.L.. (2) 
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According to [4] the 90% C.L. bound is sin < 0.024. The hmit (2) is shghtly better than the 
earher (before SNO salt phase II) hmits [5, 6, 7] sin^ ^13 < 0.031 (0.055). In [8] even stronger 3a 
upper bound is quoted: sin^ 613 < 0.041. 

The combined fit of the CHOOZ, K2K and atmospheric neutrino results (without solar ones) 
gives [3] 

sin^^ia < 0.032 (0.060), 90% (3a) C.L.. (3) 

(See also [9].) Comparison of the bounds (2) and (3) shows that effect of inclusion of the solar 
neutrino data in the global fit is weak but not negligible: the bounds improve by ~ 25%. 

Future long baseline experiments are geared towards precision measurement of 613. MINOS [10] 
and CERN to Gran-Sasso [11] alone will be able to only slightly improve the limit (2) down to 

sin^e^ < 0.015, 90% C.L.. (4) 

J-PARC accelerator experiment [12] will have substantially higher sensitivity: 

sin^ ^13 < 0.0015, 90% C.L., (5) 

but when all correlations and degeneracies of parameters are taken into account the estimated 
sensitivity reduces to sin^^ia ^ 0.004 (90% C.L.) [13]. 

The Double-CHOOZ reactor experiment [14] will be able to improve the bound down to 

sin^ ^13 = 0.008, 90% C.L. (6) 

for Am^ = 2.4 x 10^'^ cV^. It provides a clean channel of measurement of ^13 complementary to 
accelerator measurements [15], [13]. 

Detection of neutrinos from the Galactic supernova can, in principle, probe sin^ ^13 down to 
10~^ — 10~^. The problem here is large uncertainties in the original neutrino fluxes. Essentially 
one can find whether sin^ ^^13 is larger or smaller than 10^'^ Both the upper and the lower bounds 
on sin^^^i3 may be established if sin^^^i3 is in the range 10^^ — 10^^ [16, 17]. 

Ultimate sensitivity, sin^6'i3 = lO^'' — 10^^, can be achieved at the neutrino factories [18]. 

Influence of the 1-3 mixing on the detected solar neutrino fluxes has been considered in a 
number of publications before [19] - [34], [7]. The present solar neutrino data are not precise enough 
to really probe sin^ ^13 in the allowed range (2,3), but on the other hand the bounds obtained using 
the solar data only are not much weaker than those from the reactor and atmospheric neutrino 
results. Indeed, recent analysis of the solar neutrino data alone [3] gives 

sin^ ^13 < 0.038 (0.062), 90% (2a) C.L., (7) 

and very weak 3a bound. Inclusion of the latest KamLAND data improves the latter bound 
substantially [3] 

sin^^ia < 0.035 (0.075), 90% (3a)C.L.. (8) 
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(This can be compared with (3).) 

The 1-3 mixing taken at the present upper bound modifies the allowed region of the parameters 
sin^ 9 12 — Am|j obtained from the global fit. It shifts the region toward larger Am^^ and smaller 
sin^^i2, however the shift of the best fit point is within la [7]. 

It was noticed that the Earth matter regeneration effect has strong dependence on ^13, as 
P oc cos^6'i3 [35]. This can be used to determine 6*13 in future solar neutrino experiments [35, 36]. 
Measurements of the pp-neutrino fiuxes will also improve the sensitivity to sin^^i3 [21]. 

In this paper we study in details possible effects of 1-3 mixing on the solar neutrino observables. 
We analyze the present determination of 6*13 and 6*12 and evaluate the sensitivity of future solar 
neutrino experiments to ^13. 

The paper is organized as follows. We first (sec. 2) re-derive the i/g— survival probability 
including all relevant corrections and estimating accuracy of the approximations. We consider 
dependence of the probability on sin^^i3. In sec. 3 we find simple analytic expressions for the 
probability in the high and low energy limits. In sec. 4 we introduce the sin^ 612 — sin^ ^^13 plots 
and derive the equations for contours of constant values of observables in this plot. In sec. 5 
we find the allowed region of parameters sin^ $12 and sin^ ^13 and consider the bounds on these 
parameters from individual measurements. In sec. 6 the sensitivity of future, in particular, the 
low energy (LowNu) solar neutrino experiments to 1-3 mixing, is estimated. We discuss possible 
implications of future measurements of ^13 in sec 7. Our results are summarized in sec. 8. 



2 The survival probability and ^13 

The analytic formula which describes conversion of the solar neutrinos in the case of three neutrino 
mixing and mass split hierarchy, Am|i » Am|i, has been derived long time ago [31, 32]. Here 
A.mli and Am^i are the mass splits associated to 1-3 mixing and to the main channel of the solar 
neutrino conversion correspondingly. The formula was analyzed and further corrected recently [33, 
34, 35]. 

In view of future precision measurements one needs to have an accurate formula for the prob- 
ability with all relevant corrections included and errors due to approximation made estimated. 
For this reason wc rc-dcrivc the formula in a simple and adequate way which employs features 
of the Large Mixing Angle (LMA) MSW solution. We will use a very high level of adiabatic- 
ity [37, 38, 39, 40] in whole Su system. 

Due to loss of coherence between the mass eigenstates, i/i (i = 1, 2, 3), the survival probability 
can be written as 

ee- 2^ ^ei Ue ' K^) 
1=1,2,3 

where P^^"' is the probability of conversion in the Sun and p^"''^^^ is the Vi — > oscillation 

probabihty inside the Earth. Strong adiabaticity of the neutrino conversion in the Sun leads 
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immediately to 

nr = |t^e7T, (10) 

where is the ei-element of neutrino mixing matrix in matter in the production point. Indeed, 
C/^° = (I'eWmi) determines the admixture of the eigenstate, ly^i, in the neutrino state in the 
initial point. Then due to adiabaticity transforms to z/j when neutrino propagates to the 
surface of the Sun. Corrections due to the adiabaticity violation are negligible [40]: APce/Pee ~ 
7^ cos26'i2/4sin^ 9i2, where 7 is the adiabaticity parameter. For E = 10 MeV we find APee/Pee ~ 

Combining (9) and (10) we obtain 

Pee= E l^eTT^r*'- (H) 
1=1,2,3 

In the absence of the Earth matter effect (the day signal) we have 

pEarth^^^^,^2^ (12) 

where U^i are the elements of the vacuum mixing matrix. 

In what follows wc will neglect the Earth matter effect on the 1-3 mixing which is determined 
by 2EVe/ Aml^ < 0.005, (here Ve is the typical potential inside the Earth). Therefore 

Pie"''" ~ \Ue3\'. (13) 

The Earth matter effect on two other probabilities can be described by the term F^^g defined as 
the deviation of probability p^"'"*'' from the no-oscillation one: 

P2f'-*'^ = |C/e2r + F,e,. (14) 

Then the unitarity condition, Ei=i,2,3 P^"''^^ = 1 and (13) lead to 

Pir'''=\Uel\^-Freg. (15) 

Plugging (12, 13, 15) into (11) we obtain 

Pee= E \U'S'?m' + Fre,{\U:ir-m'?)- (16) 
1=1,2,3 

We will use the standard parameterization of the mixing matrix in terms of mixing angles: 

|C/el| = COS^13COS^12, | C/e2 1 = COS ^13 siu 6'i2 , | C/e3 1 = siu 6*13 , (17) 

and similar definition holds for the matrix elements in matter with substitution 9iz — ^13° and 
^12 — The angle 612 is identified with the "solar mixing angle". 
In terms of mixing angles the probability (16) can be rewritten as 

Pee = COS^ COS^ OizPad + siu^ siu^ ^13 - Freg COS^ COS 2^^^°, (18) 
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where Pad is the usual two neutrino adiabatic probabihty: 

Pad = sin' 012 + cos' cos 2^12. (19) 
For Freg = the probabihty (18) coincides with the one obtained in [32]. 

Let us find the mixing angle 9^ in matter. The Hamiltonian of system can be written in 
the flavor basis as 

H - U^zU^^Ur^H^Ul^Ulul^ + V^. (20) 

Here = diag{0, AmlJ2E, Aml^/2E), V3 = diag{V,0,0), and V = V^GppYe, where p is the 
matter density and is the electron number density fraction. 

Performing rotations U23U13 we arrive at the basis (fg, z/g) ^^3) in which the Hamiltonian becomes 

H = Ui^H'^Ul^ + UlVUi3. (21) 

If the matter effect is small we can neglect mixing of the third state (off-diagonal terms of the 
matrix Ul^VUis). This state then decouples and for the rest of the system the Hamiltonian 
becomes 

H2 = Ui2H^Ul2 + cos' ^13^2, (22) 

where V2 = diag{V,0). 

In the lowest approximation the matter effect on the 1-3 mixing can be found making an 
additional 1-3 rotation which ehminates the 1-3 elements in (21). The angle of rotation, ^^3, is 
given by 

^ Amy{2E) - Vcl2el+\ml,/{2E) sin' 6,2 ^^^^ 



where 



_ _ / E \( pYe \ /2.5 • 10-3eV'\ 

= A;;^ - I lOAfeV J I lOOiT^ ) \ ^rnl, ) ' ^^^^ 



So, approximately, — ^13613. 

The additional 1-3 rotation generates neghgible 2-3 elements ^ 



F . . Am' 

It also leads to corrections to the 1-2 block of the Hamiltonian: 



H23 = ^ sin 2^13 sin 2^12 -r^. (25) 
4 Ami 



H2 = Ui2H'^Ul2 + cos' ^13 (1 - Sin' ^13613) V2. (26) 



^Elimination of this term requires an additional 2-3 rotation by the angle ^23 = 'iEH23/ Am^i < 10 
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In the first order in 2VE/A'ml2 the rotation decouples the third state. So, the 1-3 mixing 
angle in matter is given by 

0T3^0i3 + 0[s + O{e[,e^s), (27) 

and using (23) we obtain 

sin^ « sin^ (1 + 2613) + O (sin^ O.^el^, sin^ ^13613) ■ (28) 

(See also [33].) For typical energy £■ = 10 MeV the relative matter correction is (10 - 15)%. The 
correction is negligible at low energies. 

The 1-2 mixing in matter is determined from diagonalization of (22) or (26). Neglecting very 
small correction (~ €13 sin^ Ois) in (26) one obtains for 9^ the standard 21/ formula with substitu- 
tion V cos ^13. 

Let us find dependence of Freg on the 1-3 mixing explicitly. The probability of 1/2 — > Ve 
conversion can be written in terms of the amplitudes of transitions between the mass states 
1/2^ J^i , A2, as 



2e 



AuUel + A22Ue2\ = COS^ ^13 | ^2 COS ^12 + ^22 Siu ^12 1 . (29) 



Since in the absence of matter effect Au = and ^422 = 1, in the lowest approximation in ^412 oc V, 
we obtain 

Freg = cos^ ^13 siu 2^i2-ReAi2. (30) 

The transition amplitude A12 should be found by solving the evolution equation with the Hamil- 
tonian (26). It is proportional to the potential: ^412 oc Fcos^^ia, and consequently, Freg (30) can 
be rewritten as 

Freg^COsH,sfreg + 0{f^J, (31) 

where 

freg = sin 29i2ReA[^ (32) 

is the t/g^icgcneration factor given by the 2i/ amplitude A[2 = ^12/ cos^ 6*13, and it does not 
depend on 613. 

The regeneration factor integrated with the exposure function over the zenith angle can be 
estimated as 

. £;yosin^2gi2 

^ 4Ami, ' ^^^^ 

where Vq is the potential at the surface of the Earth. Indeed, the integrated effect is mainly due 
to the oscillations in the mantle of the Earth along the trajectories which are not too close to the 
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horizontal ones. (For the core crossing trajectories the effect of the core is attenuated [41].) In 
turn, for the mantle trajectories the adiabaticity is nearly fulfilled [40] and therefore characteristics 
of oscillations (averaged probability and depth of oscillations) are determined by the potentials in 
the initial and final points of the trajectory, that is, by the potential near the surface of the Earth 
(for numerical consideration see [42]). 

Inserting the expression for F^eg into (18) we obtain 

Pee = COS^ COS^ Oi^iPad - COs' ^13 COS 29"^^ freg) + Siu^ 6^^ Siu^ ^13. (34) 

Taking the mixing angle 9"!^ in the first order in eis we find 

Pee ~ cos^ ^13(1 - 2 tan" ^laeia) {Pad - cos^ cos 29^^ Ueg) + sin^ ^13 . (35) 

Notice that the pre-factor (1 — 2 tan^ 9i'iti'i{E)) which describes the matter effect on the 1-3 mixing 
differs in the high and low energy limits. 

The expression (34) can be rewritten as 

Pee ^ cos^ cos^ 9,sP^^\Am'' , 9,2, 9is) + sin^ 9Ts sin' 9,s, (36) 

where 

P^!^ = Pad - COS^ COS 29T2'freg 

is the total Ve survival probability for 2z/ mixing which includes both the conversion in the Sun 
and the oscillations in the Earth. Pad depends on 9,3 via the effective matter potential Fcos^ ^13. 

The probability should be averaged over the distribution of the neutrino sources in the Sun. 
For the fT-component of the solar neutrino spectrum [K = pp, Be,pep, B, N,0) the averaged 
probability can be written as [40] 

Pg^fi^ = sin^ 9 + cos^ 9m{VK) cos 29 - 0.5Sk cos 29m{VK) cos 29 - cos 29m{VK) cos^ ^13/reff, (37) 

where Vk is the effective value of the matter potential in the region of production of the K neutrino 
component. Sk is an additional correction due to the integration over the production region [40], 
and in what follows we will neglect it. 

Notice there are three features of dynamics of propagation which allow one to immediately 
derive the Su survival probability in terms of the 2u probabilities: 

1) . Adiabaticity of propagation of the third neutrino; 

2) . Negligible Earth matter effect on 1-3 mixing; 

3) . Averaging of oscillations and loss of coherence of the third mass eigenstate. 
Essentially this means that the third neutrino eigenstate decouples from dynamics and propagates 
independently. Its contribution to the probability adds incoherently. Therefore the problem is 
reduced to two neutrino problem and additional factors in (34) are just projection of the fiavor 
basis onto the basis of states which contains the decouphng state us^. In (36) factors sin^ ^^5° 
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cos^ 9"^ are from the projection in the initial state, whereas the factors sin^ 9is and cos^ 9is are 
from projection in the final state. 

Up to negligible corrections the expression for the probability (36) can be rewritten in form 
Pee = (1 - sin^ 9l,fP^f{Am', 9u, hz) + sin^ hz, (38) 

where 

Sin2^i3 = sin2^i3(l + ei3), (39) 

which coincides with the form without matter corrections to ^^13 ^. This means that the matter 
effect on the 1-3 mixing in the probability can be accounted by renormalization of sin^ ^13. In- 
versely, using usual formula for the probability without matter corrections in the analysis of the 
data we determine ^13 averaged over the relevant energy interval. Then the true vacuum mixing 
angle will be smaller: 

sin2^i3 = sin2^i3(l-ei3). (40) 

Here €13 = (€13) e is the value averaged over the relevant energy range. Renormalization depends 
on the neutrino energy: 613 ~ 0.05 at high energies and 613 ~ 0.002 in the pp-neutrino range, and 
the latter can be neglected. 



3 Two limits 

As we will see, the expected sensitivity of solar neutrino studies to the 1-3 mixing will not be 
better than a{siv? ^13) ~ 0.01. According to (36) that would correspond to the change of the 
survival probability 

5P ^ -2P(7(sin2 ^13) ~ 1%. (41) 

Therefore in the present discussions we can neglect various corrections to the probability which 
are much smaller that 0.01. In particular, 

• the term sin^^i3 < 0.002 can be neglected in eq. (36); 

• the matter effect on the 1-3 mixing is small and can be taken into account as a small 
renormalization of sin^^i3 (see sec. 2); 

• we will neglect f^^g < 10~^ corrections in Freg- 

^Strictly, this is correct in the high energy or low energy hmits, where Pge^ ~ Pad depends very weakly on the 
potential. At low energies the matter corrections are negligible anyway (see eq. (24). 
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In qualitative consideration we will use simplified expressions for the survival probability ob- 
tained in the limits of high and low energies. In fact, all relevant fluxes are either in the low (pp, 
Be, N, O, pep) or in the high, E > 5 MeV (B, hep), energy limits. The expansion parameter 
equals 



at high energies and l/jy at low energies 
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1). In the high energy limit, performing expansion of (37) in rj we find the survival probability 

P^^^ ^ cos" 013 sin' 012 + \ cos 2^12 sin' 2^1277' + COS^ ^ig/.eg- (43) 

Here the first term is due to the non-oscillatory adiabatic conversion. In the second term, which 
describes the effect of averaged oscillations, dependence on ^13 is absent. In contrast, the third 
(regeneration) term has the strongest dependence on ^13 [35]. It originates from the common 
factor cos^ ^^13 and the effective potential in the 3z/ case which contains another cos' ^13 factor. 
In the lowest order in sin' ^^13 we can rewrite (43) as 

Pf^ ^ sin' 0^2 + \ cos 2^12 sin' 20^2rf + !reg - sin' ^13(2 sin' 0^2 + 3/,eJ . (44) 

Neglecting the oscillation corrections (second term) and the double suppressed term ~ sin' 0\zfreg 
we obtain 

P^^ (1 - 2 sin' ^13) sin' 0^2 + freg- (45) 



2). For the low energy neutrinos {E < 2 MeV), the matter effect is small and in the lowest 
order in r]~^ we have 

PiP K cos^ ^13(1 - 0.5 sin' 2^12) - 0.5 cos^ ^13 cos 2^12 sin' 2^1277"^ (46) 

Here the first term is the effect of averaged vacuum oscillations. The second one gives the matter 
effect correction which contains strong dependence on 6'i3. An additional second power of cos ^13 
originates from the matter potential. The regeneration effect is neglected. Numerical calculations 
confirm that the linear in rj^^ oc E approximation (46) works well up to 2 MeV. 
In the lowest order in sin' ^13 we find from (46) 

PiP 1 - 0.5 sin' 2^12 - 0.5 cos 2^12 sin' 2^i2r/-^ - (47) 



sin' 



13 



2 - sin' 2^12 - 1.5 cos 2^12 sin' 2^12?;-^] . (48) 

For very low energies (pp-neutrino range) the matter corrections can be neglected, so that 

p{i) ^ (1 _ 2 sin' ^13) (1 - 0.5 sin' 2^12). (49) 

Prom Eqs. (43 - 48) we conclude the following. 
^Notice that in [40] we use different definition of r/: r/ — > I/tj'. 
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• The dependences of the probabihty on 612 for the high energies and low energies are opposite: 
Pg^P increases whereas Pjg^ decreases with decrease of Bu- 
rn In contrast, Pj^^J and Pj^^ both decrease with increase of ^^13. So, 612 and ^^13 correlate at low 

energies and anti-correlate at high energies. For the neutrinos an increase of 613 can be 
compensated by increase of 612, whereas at low energies (for the pp neutrinos), 9i2 should 
decrease to compensate increase of ^13. 

• Both for the high and low neutrino energies the dependence of probability on Ani2i appears 
via small corrections only, and therefore, is weak. 

Apparently strong dependence on ^13 cancels in the ratio of probabilities (or the ratio of the 
corresponding observables) : 



P^^^ _ Sin^ 012 + 0.25 COS-^ ^13 cos 2^12 Sin^ 2^1277^ + COS^ Oisfreg 

W ~ 1 - 0.5 sin^ 26 



(50) 
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Notice that dependence of the ratio on ^13 appears via small corrections and it can be neglected in 
the first approximation. So, measurements of the ratio gives sin^ 9i2. Then one of the probabilities 
(in low or high limits) can be used to find sin^ ^13. 

4 sin^ O12 — sin^ Ois plots and iso-contours of observables 

There are two important features of the problem: 

• Weak dependence of the probabihty and observables on Am|i within the allowed region; 

• sin^ 6*12 — sin^ 6*13 degeneracy. 

Let us consider these features in some details. 

1). The KamLAND spectral data has shown high precision in determination of Am^i. The 
relative error 

is about 5a = 12% {3a) at the present, and with 3 kTy statistics from KamLAND this can further 
reduce down to 7% [43]. 

As follows from (43) and (48), variation of the probability in the high energy limit with Am|^ 
equals 



5PW 



^Tj'^ COS 2012 Sin^ 26I12 - COS^ Oiafreg 



5a- 



(52) 
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For E = W MeV and 5a ~ 0.1, we find 6P ~ 0.005. 



In the low energy part, variations of the probabihty with Am2i are even weaker: 

^Pee ~ 7;V~^ cos^ ^13 COS 2012 sin^ 2^i25a. (53) 
For E = 0.4 MeV we find ^Pj^) ^ q.OII^a ~ 0.001. 

2). According to the discussion in sec. 2, in the first approximation, all experiments which 
probe the high energy part of the spectrum are sensitive to the same combination of mixing angles 
6i2 and ^13. Similar statement is true for the low energy experiments but the combination of 
mixing angles is different. Therefore to determine sin^ ^12 and sin^^ia separately one needs to 
employ both types of experiments. 

In view of these two features, we will present results of our study in sin^ $12 — sin^ ^13 plane for 
fixed values of Am^i- 

Let us construct the contours of constant values of various observables in the sin^ 6*12 — sin^ 6*13 
plane. Essentially, the iso-contours for a given observable X coincide with the contours of constant 
survival probability averaged over the relevant energy range with the corresponding cross-section: 

{Ree)x = Cx. (54) 

The survival probabilities found in the high and low energy limits determine two different types 
of the iso-contours. For the high energy observables {E > 5 MeV), using Eq. (44), we find an 
analytic expression for the contours of constant X: 

. 2„ Sin^ ^12 + 0.25 COS 2^12 Sin^ 2^127/1 + (/^eg)x-Cx 

2sm 012 + S{freg)x 



where 



Y A 



m 



2 



2\ 



is the averaged over the neutrino energy value of 77^ folded with the cross-section, energy resolution 
of detector, etc.. In the allowed range of sin^ 612 (with the average value of mixing denoted by 612) 
the strongest dependence comes from the denominator of (55): 

sin^ ^13 ~ A(sin2 612 -Cx), A ^ (2 sin^ 612 + Sfreg)-\ (57) 

So, in the first approximation the contours are straight lines. This reproduces well the results of 
numerical calculations presented in sec. 5, 6. Notice that sin increases with sin 9i2- 
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For the low energy observables, Y, using (48), we obtain expressions for the iso-contours as 
. 2^ 1- 0.5 sin^ 2^12 -0.5 cos 2^^i2sin2 2^12%' -C^y 

sm 9is Ki T)— , (58) 

2-sin2 2^i2 - 1.5cos2^i2sin2 2^i2r;y^ ^ ^ 



where 



In the rough approximation 

sin^ ^13 ocl - Cy - 2 sin^ ^12(1 - sin^ Q^^) (60) 

and it has minimum at sin^ d\2 = 0.5. Other terms in (58) shift this minimum to smaller values 
of sin^ 6*12. So, in the region of interest, sin'^ 612 ~ 0.3, the lines of constant values of observables 
have negative slope: sin^ ^13 decreases with sin^ $12. 



5 Constraints on the mixing parameters from the present 
experiments 

In this section we describe the results of our numerical computations. 

We first perform the global analysis of different data samples. Some details of the analysis 
follow. The CHOOZ depends on Amg^ and in our fit we vary this parameter freely in the 
range allowed by the zenith angle SuperKamiokande atmospheric neutrino data [44]. The details 
of procedure for the solar and CHOOZ data can be found in [28] while the KamLAND analysis 
follows the method outhned in [29, 6]. We use the KamLAND data from the latest version of 
paper [45] where new source of the background has been taken into account. 

For the solar neutrinos we include the data on total rates from the radiochemical experiments 
CI [46] and Ga (SAGE, GALLEX and GNO combined) [47], the Superkamiokande day-night 
spectrum data [48], the SNO CC (charged current), ES (neutrino-electron scattering) and NC 
(neutral current) rates from the salt phase-I, the combined CC-I-ES-I-NC energy spectrum from 
the pure D2O phase [49] as well as the day- night spectrum data of CC,ES and NC events from 
the salt-phase 11 [50]. 

The neutrino fiux normalization, Jb = Fb/F§^^ , where Fb is the true fiux and F§^^ is 
the flux predicted in the SSM [51], is left as a free parameter in the fit. Essentially /b is fixed by 
the NC data from the SNO experiment. For the other fiuxes and their uncertainties we use the 
predictions from the Standard Solar Model (SSM) [51]. 

In fig. 1 we show the allowed regions of the parameters sin^ 612 — sin^ ^^13 (for two different 
values of Am^2) which follow from the global fit of all available solar neutrino data as well as data 
from reactor experiments KamLAND and CHOOZ. The contours are marginalised over /b- The 
global minimum has been obtained by letting all the parameters (including Amli ) vary freely in 
the fit. With this, the values of parameters in the global minimum come at 

Amli =8.0-10-5 eV^ sin^ ^12 = 0.31, sin^ ^13 = 0.00, = 0.84. (61) 
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Figure 1: The 90%, 95%, 99% and 99.73% CL. allowed regions of mixing angles from the global 
analysis of the solar neutrino, KamLAND and CHOOZ neutrino data. The left panel is for 
Am|]^ = 8 • 10~^ eV^ and the right panel is for Am^i = 9 ■ 10"^ eV^. The global best-fit point 
obtained by letting all the parameters (including Am2i) vary freely is marked by the crossed circle. 
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Ami, = 8.0xlO"^eV^ Am^„ = 9.0xlO"^eV^ 
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Figure 2: The 90%, 95%, 99% and 99.73% C.L. allowed regions of mixing angles from the global 
analysis of the solar neutrino data only. The left panel is for Am|^ = 8 • 10~^ eV^ and the right 
panel is for Am|]^ = 9 • 10~^ eV^. The global best-fit point obtained by letting all the parameters 
vary freely is marked by ©. Also shown by crosses are the local best-fit points obtained at fixed 
Amni . 
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The two panels in fig. 1 correspond to two fixed illustrative values of Am^^ from the current 
allowed range. The contours are plotted with respect to the global minimum (61) using the defini- 
tion of Ax^ which corresponds to 3 parameters. The upper bound, sin^ ^13 < 0.055 (99% C.L.) is 
tighter than the one obtained without including the phase-II SNO salt data (0.061). According to 
fig. 1 with increase of Am|]^ the allowed region shifts to smaller sin^ 9i2', with increase of sin^ ^13 
the value of sin^ 9i2 in minimum of practically does not change. 

We have performed the global fit of the solar neutrino data only. The best-fit values of param- 
eters are found to be 

Amj^ = 6.4 X 10"^eV^ sin^ 612 = 0.33, sin^ ^13 = 0.01, fs = 0.84. (62) 

Now the 1-3 mixing is non-zero but statistically insignificant. The solar is somewhat fiat in the 
region sin^ ^13 = - 0.03. For sin^ ^13 = 0.01 and 0.02 we obtain the xLn = 113.84 and 114.04 
correspondingly. 

Notice that we find small boron neutrino flux: fs — 0.84 in the units of flux predicted in [51]. 
Recent calculations of fluxes with new radiative opacities and the heavy element abundances give 
even smaller value fs = 0.78 [52]. 

In fig. 2 wc show the constant contours obtained from global analysis of the solar neutrino 
data only and when Am|j^ is fixed. The contours have been calculated with respect to the global 
minimum (62) and using the Ax^ values corresponding to 3 parameters. We also show the local 
minima found for fixed values of Aml^. In particular, for Am|i = 8 • 10~^ eV^ we obtain 

sin^ ^13 = 0.007 tofo7 (90% C.L.). (63) 

The best fit value of ^13 changes with Am^i very weakly, e.g., for Am2i = 9 • 10~^ eV^ we find 
sin^^i3 = 0.01. The fig. 2 shows also that at present the solar data alone has a weak sensitivity 
to sin^ 6'i3: the 90% C.L. bound equals sin^ ^^13 < 0.087. Inversely, the values of sin^ 9is allowed 
by the global fit (2) have weak impact on the global analysis of the solar neutrino data. Prom fig. 
2 we find that with increase of sin^ ^13 from to 0.05 the best fit value and 90% allowed region 
of sin^ $12 shift to larger values by A sin^ 9i2 — 0.02 only which is smaller than la error and in 
agreement with earlier estimations in [22]. 

Let us now consider bounds on mixing angles 6*12 and 6*13 from individual observables. 
1). At present only the gallium experiments can be considered as the low energy experiments. 
For the Ge-production rate, Qoe, we have 

where Qq^ is the contribution to the rate from the ith component of the solar neutrino fiux 
according to the SSM, PJg is the corresponding effective (averaged over the energy range) survival 
probability. The rate has the contributions both from the low and high energy parts of the 
spectrum, though the latter (from boron neutrinos) is much smaller. We can subtract the boron 
neutrino effect, Qoe^^, using experimental results from SNO: Qq^ = Qoe — Qat^^- Taking the 
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flux from SNO as an input, we find tliat tlie neutrino contribution equals Qq^'''^ = 4.2^07 
SNU and then = 63.9l|^ SNU for the rest of the neutrinos [47]. 

The conversion effect on Qq^ is described by the probability in low energy limit: 

Qgc = (QgT' - Qce) [cos^ ^13(1 - 0.5 sin^ 2^12) - cos^ ^13 cos 2^12 sin" 2^i2r/Ge] , (65) 

where 

nssM nB - iv;;Q'^e + V^elQ'o! + VbIQgI + VBNoQGe) ■ (66) 

At the same time we find that dependence of the total rate Qce on fs and the high energy 
conversion effect is very weak: 10% variations of fs would produce change of Qce about 0.4 SNU. 
Therefore results of numerical calculations are given in figs. 3, 4, 5, for Qce and Jb = 1. We show 
the contours of constant total rate Qce which correspond to the present central value as well as to 
la and 2cr deviations from it. The contours are well described by the analytical expression given 
in (58) with rjy^ = rj^l and Cy = Cce- 

For zero 1-3 mixing the central value Qce = 68.1 SNU corresponds to sin^ 9^ = 0.34, and the 
best fit value sin^ 6*12 = 0.33 is accepted at 0.5a level. For a given Qce, sin^ 6*13 increases with 
decrease of sin^ 612. Taking the best fit value sin^ 612 — 0.33 (see eq. (62)) and Qoe — 68.1 SNU 
we find sin^ ^13 = 0.004. 



2). The CC/NC ratio of events at SNO is given by the integrated (with cross-section) survival 
probabihty at high energies: CC/NC = {P^f')- The iso-contours are described approximately by 
Eq. (55) with Cx = Csno — CC/NC. The results of exact numerical calculations are shown in 
fig. 3. The iso-contours are constructed for the present central value of CC/NC as well as for its 
la and 2a deviations. At sin^ ^12 = 0.33 the central value of CC/NC gives sin^^i3 = 0.046. 

As seen from fig. 3, the combination of SNO CC/NC and Ga-results prefers the non-zero but 
statistically insignificant 1-3 mixing. The iso-contours of the central observed values of CC/NC 
and Qoe cross at 

sin^ = 0.017 ± 0.026 (67) 

which is about 0.7a deviation from zero. Here la is evaluated from the fig. 3 using CC/NC and 

Qoe measurements. 

The CC/NC ratio is free of the neutrino normalization factor /g. Iso-contours for the Ge- 
production rate are constructed for /s=l, i.e., for the SSM value of the fiux. The dependence 
of Qat on fs essentially disappears if one uses for the boron neutrino fiux the value measured by 
NC in SNO or by SK. In any case the influence of uncertainty related to fs on the intersection of 
iso-contours is weak. Larger uncertainty (~ 1 SNU) follows from the Be-neutrino flux contribu- 
tion. That leads to Asin^ ^13 ~ 0.01. Future Borexino measurements can reduce this uncertainty. 



3). The SuperKamiokande rate of the i/e— events. The ratio of the observed to SSM predicted 
event rates is given by 

RsK^fB[{Pee){l-r) + r], (68) 
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CC/NCvsGa 




Figure 3: The iso-contours of the CC/NC ratio at SNO and Ge-production rate in the sin^6'i2 — 
sin^ ^^13 plane for Am^i = 8 ■ 10~^ eV^. Shown are the iso-contours for the central values of 
observables and the ±lcr and ±2a deviations. We take Jb = 1- The solid (blue) contour bounds 
the 99% CL. allowed area from the global fit of all data. 
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where r = Oy^joi,^ is the ratio of the V(. — e and v^ — e cross-sections. The hnes of constant Rsk 
coincide with the hnes (Pee) = Csk — const, where according to (68) 

CsK = (69) 

1 — r 

The hnes are given by the Eq. (55) with Cx = Csk and the averaging (integration) {...)x should 
be done according to the SK experiment characteristics. 

The iso-contours of the SK rate in the sin^ 6*12 — sin^ 6*13 plane are presented in fig. 4 for 
/b = 1. They correspond to the SK central value Rsk — 0.41 as well as the la and 2a deviations. 
According to fig. 4, for sin^^ia = the ratio Rsk — 0.41 is achieved at sin^^i2 ~ 0.21 which is 
about 2a below the best global value. In turn, the best global value sin^ 612 — 0.33 and Rsk — 0.41 
lead to sin^ ^^13 = 0.17. 

With decrease of Jb the corresponding iso-contours shift to larger 6*12 and smaller 6*13. For 
arbitrary Jb one can still use the grid of contours calculated for fB — 1, but with changed 
values Rsk- Indeed, according to (69), for given values Rsk and /s one should take the contour 

Rsk = Rsk/ fs- 

The combination of SK- and Ga- data (crossing point of the contours which correspond to 
the central experimental values) selects sin^ ^^13 = 0.053 as is seen from fig. 4. For = 0.84 the 
best fit SK value, Rsk = 0.41, corresponds to the iso-contour with = 0.49. The intersec- 
tion of this contour with the iso-contour Qoe — 68.1 SNU occurs at much smaller 1-3 mixing: 
sin^ ^13 = -0.004. 



4). The Argon production rate in the Homestake experiment is determined by 

Qm = [P!:Qfr + n^QZ + n^^QSr + P^fsQl] . (70) 

According to the SSM Qat is dominated by the high energy (boron) neutrinos. Therefore the 
contours of constant rate (fig. 5) have typical slope of the high energy observables. We show these 
contours for /b = 1. 

For sin^^i3 = and the best global value of sin^ ^12 the rate Qat is about 3.15 SNU which 
is 2.5(7 above the experimental result. The combination of CI- and Ga-results prefers large 1-3 
mixing: sin^ ^13 = 0.1 which deviates from zero by ~ 2.5(7. 

For the arbitrary Jb again one needs to rescale Qat ~^ Qat/Ib- Thus, for Jb = 0.84 and 
Qat — 2.55 SNU the iso-contour with Qat ~ 3.04 SNU should be used. It intersects with iso- 
contour Qce = 68.1 SNU at sin^ ^13 = 0.032. 



Summarizing, for Amgi = 8 ■ 10^^ eV^, the combination of SNO and Ga data selects sin^ Q\j, = 
0.017, which is practically independent of Jb- For Jb = 0.84 (which follows from the global fit) 
the combination of SK and Ga results leads to sin^ ^13 = —0.004 and the combination of Gl and 
Ga data selects sin^^i3 ~ 0.032, however statistical significance of the latter is low. The average 
from these three combinations is about sin^^i3 ~ 0.01. It agrees well with results of fit of all the 
solar data for the same Am|]^ . For larger Jb the larger values of sin^ ^13 would be preferable. 
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SKvs Ga 




Figure 4: The iso-contours of Rsk and Ge-production rate in the sin^ 612 — sin^ ^^13 plane for 
Am2i = 8 ■ 10~^ eV^. Shown are the contours for the central values and ±lcr, ±2(T deviations 
from these central values. We take fs = I- The solid (blue) contour bounds the 99% CL. allowed 
area from the global fit of all data. 
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ClvsGa 




Figure 5: The iso-contours of the Ar- and Ge- production rates in the sin^ 612 — svo? ^13 plane for 
Am2]^ = 8 ■ 10~^ eV^. Shown are the contours for the central values as well as ±lo" and ±2cr 
deviations. The contours have been calculated for = 1. 
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Notice that available spectral information does not give strong restriction on sin^6'i3. In fact, 
flatness of the measured spectrum of events at SNO and SK may even favor non-zero sin^^i3. 
Indeed, at high energies an increase of sin^ ^13 implies increase of sin^ $12 and therefore flattening 
of the spectrum. 

The results on 1-3 mixing described above are in agreement with the best fit value from the 
global analysis of all solar neutrino data, sin^ ^13 = 0.01, which is reahzed, however, for smaller 
Am2i (62). Inclusion of the CHOOZ result diminishes sin^^i3 to practically zero. 

Let us consider dependence of the iso-contours on Aml^. Using (55) we find for the high energy 
part: 



^(sin^^ia) ^ 



sin^ ^13 + 



Cx - sin^ ^12 - S{frey)x{0.5 - sin^ ^13) 



Sin^ 612 + l.5{freg)x 



Sa, (71) 



where we have expressed rj in terms of sin^ ^13. 
In the low energy region according to (58) 



2 „ . riY^ cos 2^12 sin^ 29i2 



^ 4 1 - 0.5 sin^ 2^12 ^ ^ 

and the pre-factor is small: 5(sin^6'i3) ~ 0.0115a- 

In fig. 6 we show results of numerical computations of the iso-contours for Ga-rate and CC/NC 
ratio for three different values of Am^^. In agreement with our qualitative consideration a shift 
of the low energy contours is much weaker. Variations of Amji from 7 to 9 ■ 10~^ eV^ (which is 
about 25% change) gives S{sin^ Oi^) = 0.008 and 6{siia^ 612) = 0.012. Reduction of the uncertainty 
in Amli down to 14% will result in 5(sin^6'i3) = 0.004 and 5(sin^6'i2) = 0.007. This is much 
smaller than the expected sensitivity of future solar neutrino studies. The contours have been 
obtained for — 1. We have found that 10% variations of fs produce shift of the Ge-lines by 
5(sin^^i2) = 0.004, and the SNO contours do not change at all. 



Let us consider the effect of matter corrections to the 1-3 mixing. As we have established 
in sec. 2, the effect is reduced to renormahzation of sin^^i3 according to Eq. (40). True value 
of sin^ ^13 is smaller than the one extracted from the data without taking into account matter 
correction. The renormahzation is negligible for the low energy observables and it is of the order 
5% at high energies. So, the iso-contours of high energy observables should be shifted down by 
factor ~ 0.95. This, in turn leads to decrease of sin^ 6*13 in the intersection points: e.g., instead of 
sin^ ^13 = 0.017 one finds 0.016. 

The figures presented in this section allow us to estimate relevance of the non-zero 1-3 mixing 
for the solar neutrino. 

For a given CC/NC ratio an increase of sin^ ^13 from to 0.05 (~ 3a upper bound) corresponds 
to decrease the Ge-production rate by 2a or AQa^ ~ —8 SNU. Inversely, for a given Qce, the 
same increase of sin^ ^^13 leads to 1.8cr decrease of the ratio CC/NC: A CC/NC ~ —0.05. The 
same change of sin^ ^13 diminishes the SK rate by Aa: A.Rsk — —0.04. So, the la variations of 
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CC/NC vs Ga 




Figure 6: The iso-contours of CC/NC = 0.34 at SNO and Qce = 68.1 SNU in Ga experiments in 
the sin^ 612 — sin^ ^13 plane for different values of Amgi : Amji = 9 • 10~^ eV^ - the dotted lines; 
Am2i = 8 ■ 10"^ eV^ - the solid lines; ^m^i = 7 • 10~^ eV^ - the dashed lines. 
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6*13 produce about la or smaller variations of observables. As we have discussed in the beginning 
of previous section, the impact of ^13 on the global fit of the solar neutrino data is much weaker: 
3(7 increase of sin^ 613 leads to only 0.7a increase of sin^ 612. 



6 Sensitivity of the future solar neutrino experiments to 
1-3 mixing 



As we have shown, some sensitivity of the solar neutrino data to the 1-3 mixing appears essen- 
tially due to the combination of the SNO- and SK- results from the one hand side and the Gallium 
results from the other side. Future high precision measurements will have better sensitivity. To 
evaluate this sensitivity and to study possible implications of measurements we have constructed 
fine grids of the iso-contours of various observables. We will use the expected la errors, ai^, in 
determination of sin^ ^13 as a measure of the sensitivity. 2a and 3a errors can also be estimated 
from the grids of iso-contours we present. In the first approximation na error is given by nai^. 

1). In fig. 7 we show the fine grid of iso-contours for the CC/NC-ratio and the Ge- production 
rate. The SNO (phase III) will measure the NC- event rate with 6% accuracy [53] . If we assume 
a reduced error of 5% for the CC-data then this makes the CC/NC error as 8% which can be 
compared to the present error of 11%. (This corresponds to reduction of the absolute error in the 
CC/NC ratio down to 0.024 from the current value of 0.035.) 

According to fig. 7, the current Ge-production rate accuracy is AQce ~ 4 SNU. Using this 
accuracy and the expected SNO-III error, we find from fig. 7 the sensitivity {la error) to sin^ 613 
in the region sin^ ^13 ~ — 0.03: 



It is only slightly better than the present one (67). 

If wc take 5% accuracy for future possible measurements of CC/NC ratio and the Ge-production 
rate error 2.6 SNU (the later is slightly larger than the present systematic error [47], and it could 
correspond to future hypothetical large mass Gallium experiment for which statistical error is 
substantially smaller than the systematic one; this can be considered as the ultimate sensitivity 
of the Gallium experiments), then according to fig. 7 



In this case the value sin^^i3 = 0.017 would have only la significance, and in the case of very 
small ^13 the upper bound wiU be sin2^i3 < 0.024 (0.048) at 90% (Sa). 

2). The Borexino and KamLAND plan to measure signal from the ^Be neutrinos. The iso- 
contours can be well described by the low energy formula (58). Taking into account contribution 
from the non-electron (converted) neutrinos we find the iso-contours, RBm-exino, which coincide 
with iso-contours {Pee) = Cbc, where Cbc = {RBm-exino - r){l - r)~^ 



ai3(CC/NC + Ga) = 0.023. 



(73) 



ai3(CC/NC + Ga) = 0.016. 



(74) 
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CC/NC vs Ga (Reduced error) 




Figure 7: The iso-contours of the CC/NC ratio at SNO (dotted hnes) and Ge-production rate 
Qce (sohd hne) with reduced steps 0.01 and 2.6 SNU respectively, in the sin^ 6*12 — sin^ 6^13 plane. 
Ami, is held fixed at 8 ■ 10"^ eV^. 
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Borexino vs CC/NC 




Figure 8: The iso-contours of Borexino rate suppression with respect to the SSM prediction in the 
sin^ 9i2 - sin^ 6113 plane for Amgi = 8 ■ 10"^ eV^. 
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The results of calculations are shown in fig. 8. The Borexino contours have similar form to 
those for Ga-rate. 

Borexino is expected to accumulate 44000 events in 4 years of time (~ 30 events/day) for 
LMA MSW solution. This corresponds to a statistical error of only 0.5% [54]. But the accuracy 
of measurement will be dominated by the background and systematics. Furthermore, the SSM 
prediction for the ^Be fiux has 10% uncertainty which will add to the total error unless one keeps 
the normalization of ^Be flux as a free parameter to be determined from the Borexino experiment 
itself with a higher accuracy. Using fig. 8 we estimate that a 5% accuracy of measurements of the 
rate at Borexino will contribute to improvement of sensitivity of the solar studies to sin^ ^13. 

3). New low energy solar neutrino experiments will be able to measure the pp-neutrino fiux 
with high accuracy [55]. Two kinds of experiments are being discussed [56]: one - using the 
charged current reactions (LENS, MOON, SIREN) and the other - using the neutrino-electron 
scattering process (XMASS, CLEAN, HERON, MUNU, GENIUS). 

For the pp-neutrinos we can use the iso-contour equation for the low energies as in (58) with 
substitution Cy = Cpp = {Pee)pp and r]~p = {EV/ Aml^)pp. Practically, r)~p is negligible. 

In fig. 9 we show the contours of constant rate of the u — e scattering events due to the 
pp-neutrinos. The neutral current contribution is also included. According to fig. 9, future mea- 
surements of the CC/NC ratio with 5% accuracy and the pp-rate with 2% accuracy will have a 
sensitivity {la error) 



Consequently, the value sin^ ^13 = 0.034 will be established with 2.3cr significance, and in the case 
of very small ^13 the upper bounds sin^ ^13 < 0.022 (0.045) at 90% {3a) can be achieved. 

4). Future Megaton scale water Cherenkov detectors HyperKamiokande [57] and UNO [58] 
will be able to measure the ratio R^e with accuracy AR^e ~ 0.003. We construct the fine grid of 
the Ri,e iso-contours for /b = 1 (fig. 10). According to fig. 10 these measurements of R,ye and the 
pp-rate (2% accuracy) would have a sensitivity {la error) 



However, the problem here is the poor knowledge of Jb- It is difficult to expect that accuracy of 
the theoretical predictions will be better than 10%. Comparable accuracy will be achieved in the 
direct measurements of the neutral currents. The global fit of the data which include the spectral 
information and the regeneration effect with free J'b may give better accuracy. 

So, we conclude that the future solar neutrino studies may reach a sensitivity sin^ 6*13 ~ 
0.015 — 0.020 {la), as compared with the present sensitivity 0.05 - 0.06. It seems that even 
combined fit of several future precision measurements will not go down to 0.01 which may be 
considered as an ultimate sensitivity of solar neutrino measurements of sin^^i3. 

On the basis of studies performed in sees. 5 and 6 we can identify signatures of the non-zero 
sin'^ ^^13. In general, they show up as a mismatch of the low and high energy measurements. The 
signatures include 



(713 (CC/NC +pp) = 0.015. 



(75) 



(713 (CC/NC +pp) = 0.011. 



(76) 
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isoppvs CC/NC 




Figure 9: The iso-contours of the CC/NC ratio (dashed hnes) and ve scattering rate of pp-neutrinos 
in the sin^ 612 — sin^ ^^13 plane for Am2i = 8 ■ 10^^ eV^. 
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pp vs HyperK 




Figure 10: The iso-contours of the neutrino electron scattering rate for the pp-neutrinos (solid 
line) and high energy neutrinos (dotted line) in the sin^ 9^ — sin^ ^13 plane. The iso-rates for 
the high energy neutrinos are drawn with an step 0.003 corresponding to the expected statistical 
error of a Megaton water detector like HyperKamiokande. Am^i is held fixed at 8 ■ 10~^ eV^ and 

fB = l. 
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• Small value of the Ge production rate Qce- For a given ratio CC/NC one can find the critical 
value of QGe(CC/NC), so that inequality 

Qoe < Q^e(CC/NC) (77) 

will testify for the non-zero value of sin^^i3. Prom the fig. 7 we find (5^g(0.30) = 76 SNU, 
g^^(0.32) = 74 SNU, Qhei^M) = 72 SNU, Q^e(0-36) = 70 SNU, etc.. The larger the 
difference in (77), the larger sin^^i3. In fact, the present data satisfy the inequality (77), 
however the statistical significance is low. Apparently if CC/NC and Qce are measured with 
high accuracy, the presence of non-zero sin^ ^13 can be estabhshed. 

The inequality (77) can be inverted, so that small values of CC/NC for a given Qce, 

CC/NC < CC/NC'=(QGe), (78) 

testify for non-zero ^13. 

• Small Borexino rate: 

Rbov < i?L.(CC/NC). (79) 

The critical Borexino rates for different values of CC/NC equal: i?^Q^(0.280) = 0.71, 
i?^„^(0.305) = 0.69, R%^{0.S25) = 0.67, i?^„^(0.340) = 0.66, R%^{0.36) = 0.645. 

• Small ue event rate produced by the pp-neutrinos: 

Rpp < R;^{CC/NC). (80) 

According to fig. 10 the critical values i?^p(CC/NC) equal 0.75 (0.294), 0.73(0.318), 0.71(0.343), 
0.69(0.373). 

• Weak spectrum distortion of the boron neutrinos. The spectrum distortion can be char- 
acterized, e.g., by the ratio of suppressions at low and hight energies vr = Rue{E < 
7 MeV)/R^^{E > 7 MeV). There is certain critical value of the distortion, r|j, which 
depends on CC/NC, so that weaker distortion will testify for non-zero 1-3 mixing. Indeed, 
for fixed CC/NC with increase of sin^^is, the 1-2 mixing, sin^^^i2, should increase (45) and 
therefore, the up turn of spectrum at low energies should be weaker. 



7 Implications of the solar neutrino measurements of ^13 

Let us first compare sensitivities of the future solar neutrino experiments and future reactor and 
accelerator experiments given in eqns. (4, 5, 6). With the expected accuracy of SNO-III and 2% 

accuracy of the pp-neutrino flux measurements the solar neutrino experiments can reach sensitivity 
comparable to the one of the forthcoming MINOS and CNGS experiments. This sensitivity is at 
the level of the present bound from the global flt of all oscillation data and it is certainly lower 
than the sensitivity of future measurements. 
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The solar neutrino sensitivity would be comparable to that of J-PARC and Double-CHOOZ if 
the pp-neutrino flux is measured with 1% accuracy and CC/NC ratio - with 3% accuracy which 
looks rather unrealistic now. 

Notice that the survival probability Pee does not depend on the CP phase S, and therefore the 
problem of degeneracy with 6 [59] does not arise. So, the solar neutrinos studies can provide a 
clean channel for determination of ^13 like the reactor neutrinos. Therefore, in spite of the lower 
sensitivity the solar neutrino data will contribute somehow to resolution of degeneracy of these 
parameters in the future global fits. 

It seems that the reactor and accelerator experiments will obtain results much before future 
high statistics solar neutrino experiments will start to operate. In this case possible implications 
of the solar neutrino measurements can include 

1) . Precise determination of sin^ Ou- resolution of the ^12 — ^13 degeneracy. 

2) . Comparison of results on ^13 (as well as on O^) from the solar neutrinos and reac- 
tor/accelerator experiments. The outcome could be further confirmation of the reactor/accelerator 
results by solar neutrinos, or establishing certain discrepancy which may imply new physics beyond 
the LMA solution of the solar neutrino problem. 

Let us consider these two points in some details. 

1) . Resolving degeneracy of angles. Due to the strong sin^ ^^12 — sin^ ^13 degeneracy unknown 
value of 9i3 leads to uncertainty in the determination of sin^ 612 and vice-versa. Thus, in the 
CC/NC measurements, inclusion of the 1-3 mixing produces an increase of sin^ 612 by 0.04 — 0.05 
(that is, by 20%) if sin^ ^13 increases from to 0.06 (see fig. 7). In the pp-neutrino fiux mea- 
surements sin^^i2 decreases by 0.06 — 0.08, when sin^^i3 increases from to 0.06 (see fig. 9). 
So, precise determination of sin^ ^12 (sin^^ia) from the solar neutrino data is not possible unless 
sin^ 9i3 (sin^ 612) is measured or strongly restricted. As follows from our analysis a combination 
of high statistics measurement at high and low energies can remove the degeneracy. 

2) . Confronting the solar and reactor/accelerator results. If non-zero sin^ ^13 is measured, our 
analysis will allow to understand relevance of the 1-3 mixing to the precision solar neutrino studies. 
For fixed Amgi, the 1-3 mixing will change determination of 612. 

Notice that dependence of sensitivity of the solar neutrino studies to 6*13 on the mass Am|i is 
very weak (via matter corrections to 1-3 mixing) in contrast to laboratory measurements. So, for 
small Am^3 the former will have some advantage. 

If 613 effect is not found in future laboratory experiments and the upper boimd at the level 
sin^ ^13 < 0.01 is established, the effect of 1-3 mixing will be irrelevant even for measurements 
with the 1% accuracy. The uncertainty of future measurements of sin^ 9^2 due to unknown ^^13 will 
be reduced down to 0.01. Furthermore, the precise measurements of sin^ 612 by CC/NC will have 
the smallest uncertainty due to larger slope of the iso-lines. 

The analysis performed in this paper will allow us to understand new physics if inconsistencies 
in results of measurements will show up. Implications will depend on character of results, and 
in the positive case, on particular value of sin^^is. Comparing the laboratory bound with the 
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results of solar neutrino studies one can realize new physics effects which may show up in the solar 
neutrinos only and can not be found otherwise. Let us consider some examples. 

Suppose that in the standard analysis precision solar neutrino measurements will give 
a nonzero value of sin^^is which is larger than the upper bound from the reactor/accelerator 
experiments. Such a situation may indicate some additional suppression of the solar neutrino 
flux at low energies. Indeed, this would shift the low iso-contours to the left - down. It would be 
equivalent of taking the low iso-contours in our plots {e.g., in fig. 3) with larger survival probability. 
The intersection with high energy iso-contours will shift down and therefore the inferred value of 
sin^ ^13 will be reduced. This situation can be reproduced in the presence of sterile neutrino with 
Am^ < 10^^ eV^ [60]. In this case the extracted value of sin^ ^^12 reduces too. 

Agreement between the solar and laboratory measurements can be also achieved by an addi- 
tional suppression of the high energy (boron) flux due to flavor conversion, so that CC/NC lines 
move to the right - down. In this case also the intersection of the high and low iso-contours 
moves down and the inferred value of sin^ $12 increases. Such an effect can be produced by certain 
non-standard neutrino interactions. 

The value of sin^ 9is extracted from the solar neutrino studies can be diminished if both the 
low and high energy fluxes are reduced by some additional mechanism. In fact, this imitates the 
effect of the 1-3 mixing and can be achieved by additional non-standard interactions. 

Let us consider the opposite situation: sin^ ^^13 is discovered in the reactor/accelerator experi- 
ments, but the solar neutrino data prefer zero or smaller value of sin^^ia. This can be explained 
if the survival probability is larger than the one predicted for given values of 6*12 and 6*13. Again 
some scenarios with modifled dynamics of conversion due to new interactions of neutrinos can 
reproduce such a situation. 

Similarly discrepancy can be realized comparing values of the 1-2 mixing extracted from the 
solar neutrino data and reactor/accelerator experiments. Indeed, the sin^^i2 can be measured 
with 3% accuracy at la (1 d.o.f.) in the reactor experiments with basehne 50 - 70 km [61, 62, 
43] provided that the 1-3 mixing is determined (restricted) accurately enough from independent 
measurements. Since a non-zero ^^13 drives the 1-2 mixing to smaller values in such an experiment, 
combination with high energy solar neutrino experiments for which non-zero ^13 drives sin^ 9i2 to 
higher values can resolve the ^^13 — 612 degeneracy. This is similar to combination of low energy 
and high energy solar neutrino experiments as discussed in the present article. In this case one can 
compare the (sin^ $12 — sin^ ^13) allowed regions from the solar neutrino studies and from reactor 
experiments. 

8 Conclusion 

1). We have derived formula for the survival probability in the three neutrino context (35) which 
includes all corrections relevant for the future precision of measurements. We estimated the accu- 
racy of approximation made. 
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2) . In the three generation context in the first approximation (neglecting the matter effect on 613) 
the solar neutrino probabilities are functions of the mass split Amli and the mixing angles: $12 
and ^13. The split Aml^ is already determined with a high precision by current data. Further- 
more, the survival probability has a weak dependence on Am|j in the allowed region. In view of 
this we have constructed various plots in sin^ 612 — sin^ ^13 plane which clearly demonstrate the 
degeneracy of these mixing angles and also show how the uncertainty in one is going to affect the 
precision determination of the other. Synergy between the high and low energy solar neutrino flux 
measurement provides a better constraint on ^13 and allows us to resolve the degeneracy of these 
mixing angles. 

3) . Among the existing solar neutrino experiments only the Ga-experiment is sensitive to the low 
energy part of the spectrum and in particular to the pp-neutrinos. Therefore we construct the 
iso-rates for Ga experiments in the sin^ 9i2 — sin^ ^13 plane and investigate how the combinations 
of these with the iso-ratcs of CI, SK and CC/NC ratio in SNO constrain ^13. 

At present the sensitivity of the solar neutrino experiments to ^^13 is rather weak and no 
indications of non-zero value has been obtained. Our analysis of all solar neutrino data leads to 
statistically insignificant non-zero best fit value: sin ^13 = 0.01. Similar value, sin^^ia ~ 0.007, 
follows from the fit of the solar neutrino data for Am^ = 8 • 10"^ eV^ selected by the KamLAND 
experiment. The combination of the CC/NC ratio at SNO and Qce rate, which are essentially 
independent of the solar neutrino fluxes, gives sin^ ^^13 = 0.017 ± 0.026 (lo"). 

The global fit of all existing data (including reactor, accelerator and atmospheric neutrino re- 
sults) leads to sin^ 613 = 0.00 as the best fit point and to the upper bounds sin^ 613 = 0.041 (0.061) 
at 90% (3(7) from the two parameter plots in sin^ 612 — sin^ 9is plane. 

4) . The precision of ^^13 determination will improve with future accurate measurements of the 
pp-neutrino flux, the CC/NC ratio and the energy spectrum of events at E > 4 MeV. We find 
that measurements with 5% error in the CC/NC ratio and 2% error in the pp-neutrino flux in ue 
scattering will have a sensitivity {la error) sin^^is ~ 0.015 (Icr). The value sin^^is ~ 0.01 (la) 
looks like the ultimate sensitivity which could be reached in the solar neutrino experiments of the 
next generation. 

Since the probability involved is Pee the problem of degeneracy with CP phase 6 does not 
appear and solar neutrinos can provide a clean channel for 613 measurements like the reactor neu- 
trinos. 

5) . It is likely that new accelerator and reactor results will be obtained first. Comparison of 
the results from solar neutrinos and accelerator/reactor experiments may confirm each other and 
further improve determination of sin^ ^13 or reveal some discrepancy which will indicate a new 
physics beyond usual interpretation of the solar neutrino anomaly. 

Determination of sin^ ^13 will eliminate or reduce the 612 — ^13 degeneracy thus improving 
measurements of 1-2 mixing. Inversely, the better determination of sin^ 612 can help in improving 
the precision of the sin^ ^13 measurements in the solar neutrino studies. 

Comparison of the results of sin^ ^13 (as well as sin^ 6*12) measurements in solar neutrinos and 
in reactor and accelerator experiments may uncover new interesting physics beyond LMA MSW 
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picture which can not be found otherwise. 
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